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Abstract 

The dynamics of an unbiased spin-boson model with Lorentzian spectral density is investigated 
theoretically in terms of the perturbation theory based on a unitary transformation. The non- 
equilibrium correlation function P{t) and susceptibility x"i^) ^-re calculated for both the off- 
resonance case A < 0.5n and the on-resonance case A ~ il. The approach is checked by the 
Shiba's relation and the sum rule. Besides, the coherent-incoherent transition point Oc can be 
determined, which has not been demonstrated for the structured bath by previous authors up to 
our knowledge. 
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I. INTRODUCTION 



Quantum computation has shown a lot of advantages in performing certain tasks 
As the basis of quantum computer, quantum bit (qubit) is one of the most attractive research 
topics today. Although, qubits have been realized in microscopic systems many years ago0, 



6| , it is difficult to implement the desired large number of interacting qubits which would 
be of practical value for computation 7|]. Therefore, macroscopic qubit systems, especially 
for the solid state circuits systems, have aroused a lot of interests recently, not only for its 
potential in realizing valuable quantum computer, but also for its theoretical importance in 
understanding the boundary between classical and quantum physics. Actually, about twenty 
years ago, macroscopic quantum phenomena has already been proposed in A.J.Leggett's 



pioneering work 



field of scalable qubits, especiall 
phase and fiux qubits 



Through several years' efforts, rapid progress has been made in the 
or the superconducting qubits including charge 



13,04 
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16|. 



>|9|,0,ll,[l2|,H,[l6|, 



Although Rabi oscillations have been observed in those qubit systemsj; 
the decoherence is still the biggest obstacle on the avenue toward realizing a quantum 
computer today. For a quantum system, any irreversible interaction with the environ- 
ment will completely destroy the quantum coherence. This raises the problem of exploring 



the dynamics o 
example 



qubit system with dissipative environment. Take fiux qubit systems for 



201 1 : during measurement, the qubit is coupled with the detector, which it- 



self is coupled to the outside environment. Via the coupling which extracts information from 
qubit system, the noise is also transmitted to the qubit system from SQUID's noncoherent 
environment. In this case, the environment affects significantly the qubit system and has to 
be taken into account. In order to explore the effect of environment on the qubit system, we 
can introduce an ohmic boson bath as a description of the non-coherent environment. And 
it has been proved that this kind of qubit system can be mapped to spin-boson model and 
the only difference is the spectral density 18, Q, [2^, so we only have to concentrate on the 



spin-boson model with different structured baths. 

In this paper, we study the dynamics of a flux qubit measured by a dc-SQUID. After 
mapping the system to spin-boson model, the corresponding spectral density is of Lorentzian 
form: 

Ji^) = 77r,. -^^T-77^ZT^-?^.- (1) 



(fi2 _ cu2)2 + (27rrcjr])2 
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We study this spin-boson model in terms of the perturbation treatment based on a unitary 



transformation which was proposed by H. Zheng[2l|. This perturbation method can lead 
to the analytical results for the non-equilibrium correlation function and the susceptibility. 
Actually, many methods have been used to study this kind of systems, however, a spectral 
density of the type Eq|T] poses a challenge to most of these approaches 241] . By comparison 
with those approaches, our perturbation method demonstrates some good features: it works 
well for a wide parameter range. And coherence-decoherence transition point ac can be 
calculated. 

This paper is organized as follows: in Sec. [Tll we present the model and give an alternative 
view of this model. In Sec. Illlt we introduce the unitary transformation briefly. In Sec. IIVI 
we calculate the non-equilibrium correlation function compare the results with that of 
other methods. In Sec. |V] the Green's function and susceptibility are calculated, and Shiba's 
relation are validated. In Sec. IVll the coherence-decoherence transition point ac is studied 
in detail. 



II. QUBIT-ENVIRONMENT INTERACTING HAMILTONIAN 

In a flux qubit system, the qubit is the two macroscopically distinct quantum states 
representing clockwise and anticlockwise rotating supercurrents. And information in the 



qubit is detected by the outside circuits including a dc-SQUID[l5|, ll6|. In this case, the qubit 
is entangled with the detecting fleld, which is itself coupled with the outside noncoherent 
environment. The qubit can be characterized by a pseudospin-1/2 operator as usual 
(unbiased condition), the detecting equipment, which is actually a LC resonant circuit 
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161], can be described by a harmonic oscillator with a characteristic frequency and the 
outside environment can be described by a set of harmonic oscillators. Therefore, the qubit- 
meter-environment interaction hamiltonian can be written as(/i =1): 

A P2 



2mfc 2 mfcCUfc 



+ E 



(2) 

where A represents the frequency of tunneling between the two states of the qubit, $7 is the 
frequency of meter, tu^'s are the frequencies of the oscillators which represent the outside 
environment [k = 1, 2, 3, ■ ■ ■), qa^ is the displacement of qubit caused by the interaction with 



meter, which also has a displacement of CkX/ {mkOJk) caused by the interaction with outside 
environment. Here, the coupling between qubit and meter is assumed to be linear, and the 
same assumption is applied to the coupling between meter and outside environment j^] • 
According to the second quantization process, the Hamiltonian H can be written as: 



+ (i?t + B) + V Kkihl + hk)] + (fit + Bf V ^, (3) 

where B (or B"^) and hk (or h\) are the annihilation (or creation) operators of harmonic 
oscillators with frequencies and cj^'s, respectively, g and Kk are the coupling constants. The 
coupling to the environment is fully defined by the spectral density,which is usually taken 
to be of ohmic form to model the dissipative environment: J{uj) = ~ ^k) = 

ruj6{u!c — uj). 



The Hamiltonian can be mapped to spin boson model 18|, |27|] : 

A 



H = -^(rx + ^ujkblbk + ^(r,^gk{bl + bk), (4) 

k k 

where the spin dynamics depends only on the Lorentzian structured spectral density which 
is approximately given by: 

T/ \ 2 c-/ \ 2aLUil^9(iUc — id) 

where a = STg'^/Q'^. Since the cutting frequency Uc is always much larger than Q, the 
spectral density J{uj) can be reasonably taken as Eq.([T]). When the characteristic frequency 
Q is higher than others, say Q > 2 A, the Lorentzian structured spectral density is nearly 
the same as the Ohmic spectral density which has been extensively studied. This case will 
be called as off-resonance. For the on-resonance case ~ A, the physical properties of the 
coupling system with Lorentzian structured spectral density may be quite different from 
those of the Ohmic bath. 

III. UNITARY TRANSFORMATION 

n 

Now we apply an unitary transformation to the Hamiltonian (jl])[2i]:if' = 

exp{S)H exp{—S), where 
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If ^jfc = 1, the generator S reduces to the usual polaron transformation. After the unitary 
transformation, the Hamiltonian can be decomposed into three parts: 



H' ^ Hq + H[ + H2, 



where 



k 



(7) 



(8) 



k 



9k 



^k 



(9) 




Y.-TUbl-bk) 
-Uhi - bk) 

^k 



-V 



Y.-^M-bk)]. (10) 



Obviously, can be solved exactly since the spin and bosons arc decoupled in H^. 
The eigenstates of Hq can be expressed as a direct product: Is)!!?^*;}), where \s) is the 
eigenstates of a^, which can be or |s2), and Ifrifc}) is the eigenstates of phonons, which 
means that there are phonons for mode k. Therefore, the ground state of Hq is given 
by : \go) = I si)|{Ofc}), where | Si) is the lower eigenstate of spin and 1(0^}) stands for the 
vacuum state for phonons. Also the lowest excited states can be denoted as | S2)|{0fc}) and 
I ■^i)!!-'-*;}) where |{lfe}) is the number state with = 1 but nk' = for all k' 7^ k. 

Since H[ and H2 will be treated as perturbation, they should be as small as possible. In 
order to minimize H[ and H2, we let H[\gQ) — and (5'o|-f^2l5'o) — 0- Then, the parameters 
77 and ^fc's are determined as. 



77 = exp 



E Sk c2 



k 



2a; 



(11) 



(12) 
1 if the 



Note that < < 1 measures the intensity of the spin-boson coupling: ^j. 
boson frequency k is larger than the renormalized tunnelling rjA; but <^ 1 for cuk <^ 77A. 
Since the transformation generated by S is a displacement, physically, one can see that 
high-frequency bosons {cuk > JyA) follow the tunnelling particle adiabatically because the 



5 



displacement is Qk^k/^k ~ Qk/^k- However, bosons of low- frequency modes ujk < 77 A in 
general are not always in equilibrium with the tunnelling particle, and hence the particle 
moves in a retarded potential arising from the low-frequency modes. When the non-adiabatic 

effect dominates, ojk <^ ^^A, the displacement <^ 1. 

The elements of Hamiltonian matrix of ground and lowest excited states can be written 

as: 

1^0) 1^0) l^fc) IV'fe') 



{i>k'\ 





2 " 



2 





Vk 



-t+ook 





Vk' 




(13) 





Vk 

Vk' 

where H2 has been dropped, \iIjq) and {ipk) (k=l,2,---) are the lowest exited states which 
represent | S2)|{0}) and | respectively. I4 = r]Agk^k/^^k = QkV^/i^k + '7A). 

Since \go) has already been diagonalized, we can now diagonalize the lowest excited states 
of H' as follows: 



(14) 



-|- terms with higher excited states . 

The transformation is given by: 

\E) = x(E)|^o) + 5^2/fc(i^)|^fc), 

k 

iV'o) = \s,)m) = Y,<E)\E), 

E 

\^k) = |si)|{u}) = ^y,(£;)|£;). 



where 



x{E) = 

yk{E) = 



Vk' 



Vk 



E + ^-cuk 



xiE), 



(15) 
(16) 
(17) 

(18) 
(19) 



where E's are the diagonalized excitation energy and they are the solutions of eigenvalue 
equation: 



Vk' 



E+^-cok 



0. 



(20) 
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IV. THE NON-EQUILIBRIUM CORRELATION FUNCTION 



The non-equilibrium correlation function P(t) is defined as P(t) = {(f>(t)\az\4>(t)), and 
\(f){t)) is the wave function in the Schrodinger picture: 

|0(t)) = 6-^^*10(0)). (21) 

we choose the initial state as \+)\b, +) where |+) is the eigenstate of cr^ = +1 and \b, +) is 
the state of bosons adjusted to the state of cr^ = +1. Because of the unitary transformation, 
the non-equilibrium correlation function P{t) can be written as: 

Pit) = ({0}|(+|e^^'V,e-^^'*|+)l{0}). (22) 
By using Eqs 11111201 we can get: 



E 



1 r f dE'e-'^'' 



Am \ J E' - TjA - R{(jj) + i-i{uj) 
dE'e'^'* 



E'-7]A- R{uj) - i-f{uj) ) ' ^^^^ 
where a change of the variable E' = E+r]A/2 is made and residue theory has been used. The 
contour of the first integrand in Eql23] is composed of a straight line which is infinitesimally 
close to the real axis from above and a semicircle above the real axis with infinite radius. 
And the contour of the second one is composed of a straight line which is infinitesimally 
close to the real axis from below and a semicircle under the real axis with infinite radius. 
R{E') and ±7(£") in Eqj23] are the real and imaginary parts of /{E'^iQ'^ — ujk) and 

they can be written as: 

7(lu) = Ti{r]Afj{uj)/{uo + 7]Af. (25) 
The integral in Eq. (l23|) can proceed by calculating the residue of integrand and the result 



is: 



P{t) = e-^^cos{ujQt), (26) 
where ojq is the solution of equation: 

uj-T]A- R{uj) = 0, (27) 



and 7 = 'y{r]A) = 7rJ{r]A)/A where we have apphed the second order approximation [28 1 . 
Once the parameters a, T and Q of the system are given, the renormahzed frequency of the 
tunnelhng between two states can be derived from Eq. (l27|) . 

For checking our approach, we calculate the renormalized tunnelling frequency c^o ioT both 
the off- and on-resonance cases and make comparison with results of the other methods. The 
main plot in Fig{T] describes the off- resonance case. We can see that the tunnelling frequency 
ujo decreases as the coupling strength a increases. This is the similar behavior as that of 
Ohmic case because the Lorentzian structured spectral density becomes the ordinary Ohmic 
one for A/Q ^ 1. But it is quite different in the on-resonant case with A ^ Q (the inset 
of Fig.2). We can see that as the coupling strength a increases, the tunnelling 



requency ujq 



increases too. It is said that the coupling enhances the tunnelling frequency 

Another way for calculating P{t), which may be more precise, is to do the investigation 
in Eq.(22) directly with the help of Kramers-Kronig relation: 

P(f) = - r duj-i{u)cos{ut) 

Trio [uJ-vA-RiiuW + ^'iujy ^ ' 

Fig 12] shows our calculation, which is compared with the result of quasiadiabatic propagator 
path-integral (QUAPI) method. Note that it is in the on-resonance case and one can see 
that P{t) shows a double-frequency oscillation. 

Figj3] shows the behavior of P{t) for different ratio A/Vt. It is a single-frequency oscil- 
lation in off- resonance case [A/Vt ^ 0.5). However, in the on-resonance region (A ~ f2), 
P{t) has two characteristic frequencies. FigJH shows the dynamics of non-equilibrium corre- 
lation function P{t) for the on-resonance case for different coupling. Its behavior changes 
significantly as a increases and becomes an over-damping curve as a becomes large enough. 

The Fourier transformation of P{t) is given by: 

^^""^ " 2^[u-r,A-R{l)f + ^'^{ujy ^^^^ 
and is shown for off- and on-resonance case in FigJSl P{t) has two characteristic frequencies 
uji and uj2i where uji corresponds to the tunnelling A and uj2 to the detecting frequency 
Vt. Since qubit is entangled with the detecting system, it is reasonable that the dynamic of 
qubit shows some feature of the detector. The results also show that: For off-resonance case 
A < 0.5i7, we have uJi < uj2 and Ui decreases as the coupling a increases. When A ~ we 
have oji > uj2 and ui increases as the coupling a increases. All these results are consistent 
with previous authors. 
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V. THE SUSCEPTIBILITY AND SHIBA'S RELATION 



The retarded Green's function is defined as: 



(30) 



where means the average with thermodynamic probabihty exp(— and [A, -B] is 

the commutator AB — BA. The Fourier transformation of G{t) is denoted as G{uj) which 
satisfies an infinite chain of equation of motion. We make the cutoff approximation for the 



equation chain at the second order of Qk and the solution for T = is 

1 

The susceptibihty xi}^) = ~G{uj), and its imaginary part is: 

7(co')6'(ci;) 



2l| 



(31) 



[uj-riA-R{uj)]^ + -f^{uj) 
^i-iu)ei-Lu) 



[LU + r]A + R{-iu)]^ + ^^{-uj) 
Define function ^(a;) as: S{uj) = 'x"{ijj)/uj with its hmit at a; ^ 0: 

hm Siu}) = - — — r^TTT- 



(32) 



(33) 



Besides, the static susceptibihty can be obtained from the imaginary part by the Kramers- 
Kronig relation: 

1 rx"iu;) 



X'iu; = 0) 

Figl6] shows that the Shiba's relation 29 



-duj. 



30, 



31 



32| 



TT 



lim^(u;) = -a[xV = 0)f 
is exactly satisfied. This fact is also a check for our approach. 



(34) 



(35) 



VI. THE COHERENT-INCOHERENT TRANSITION 



For the coherent oscillation, S{uj) has a double peak structure symmetrical with respect 
to a; = 0. However, as soon as the system becomes incoherent, S{uj) would have only a 
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quasi-elastic peak at around uj = . Therefore, the coherent- incoherent transition point 
Oc can be determined by investigating the behavior of ^(ci;). The off-resonance case where 
A/Q ^ 0.5 and the on- resonance case where A ~ ^7 will be treated separately. From FigJTl 
we can see that the coherent-incoherent transition point is at Oc = 0.49996 when A/Q = 0.5, 
r = 0.02, but it is at ac = 0.156 when A/Q = 1.1, F = 0.02. It has been validated that the 
sum rule, P(t = 0) = 1 by the integration of Eq.(27), is always satisfied when a < ac- 

A. off- resonance case: A/Q, ^0.5 

In the off-resonance case where A/Q = 0.5 and F = 0.02, according to FiglTJ we have 
ttc = 0.49997. When a goes to ac, the inherent frequency ui, which corresponds to A, 
becomes the dominant frequency of P{t) and goes to 0. In addition, from Eql28] we have: 
P(0) ~ 0.999999 for a = ac, where sum rule is still satisfied. Further more, according to 
FiglHl we can also find the dynamics of non-equilibrium correlation function P{t) is very 
like a over damping curve at the transition point a = ac- All these results show that 
coherent-incoherent transition occurs when a = 0.49996 for A/Q = 0.5 and F = 0.02. 

B. on-resonance case: A ~ 

In the on-resonance case with A/Q = 1.1 and F = 0.02, Fig|7| shows that ac = 0.156. 
When a increases to ac, the inherent frequency uj2, which corresponds to Q, becomes the 
dominant frequency of the non-equilibrium correlation function, and goes to 0. In addition, 
the integration Eq. (1251) leads to -P(O) ~ 0.999999 for a = ac where sum rule is still satisfied. 
Furthermore, Figl9] shows that the dynamics of non-equilibrium correlation function P{t) 
looks like a over damping curve. These results indicate that coherent-incoherent transition 
is at a = 0.156 for A/Q = 1.1 and F = 0.02. 

C. coherent-incoherent transition point Uc 

By calculating S{u!), the coherent-incoherent transition point ac can be determined as 
shown in FigJTOl ac = 1 /2 at the scaling limit A/Q -C 1 , which is the same as was predicted 
by previous authors in case of ohmic bath. As the system deviates from the scaling limit, 
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the coherent-incoherent transition point Oc is always less than 1/2, which is different from 
the ohmic case where ac is always larger than 1/2 for finite A. 21| 

In the off-resonance case, decreases smoothly as A/fi increases, and the descending 
becomes faster with increasing F. Whereas, in the on-resonance case, there is a suddenly 
drop of to a much smaller value at some particular values of A^/fi and r^. Our explanation 
is that as A approaches fi, the qubit-oscillator system becomes more and more on-resonant 
and the outside incoherent factor will become more and more easier to transfer into the qubit 
system under the help of oscillator fi. So, when A/fi and F reaches a particular on-resonant 
point Ac/i^ and Fc, the qubit system will show a sudden fall of coherence. FigJTTl shows that 
Fc decreases linearly as A^/fi increases. Therefore, in order to maintain coherent dynamics 
(large Oc), we have to limit F and A/fi in a particular range in the on-resonance case. 



VII. CONCLUSION 



The dynamics of an unbiased spin-boson model with Lorentzian spectral density is inves- 
tigated through a perturbation method based on a unitary transformation. An alternative 
view of the system is a two state system coupled to a single harmonic oscillator with fre- 
quency f2, the latter being weakly coupled to an Ohmic bath. By comparing with others, our 
approach shows some advantages: it works well for both the off-resonance case A ^ 0.5^ 
and the on-resonance case A ~ fi, and the coupling constant a may be as large as the 
coherence- incoherence transition point a^- We calculate the non-equilibrium correlation 
function Pit) and the susceptibility with the Shiba's relation exactly satisfied. Be- 

sides, the coherent-incoherent transition point ac can be determined, which has not been 
demonstrated for the structured bath by previous authors up to our knowledge. 
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FIG. 1: The renormalized tunneling frequency cuq verses the couphng strength a. Main plot: 
A/ri=0.1, r=0.02; inset: A/il=l.l, r=0.06. The results of the other methods are from 
F.K.Wilhem's articlej2i]. 
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FIG. 2: The dynamics of non-equilibrium correlation function P{t). The parameters are: a = 
0.004, A = r = 0.014. The result of quasiadiabatic propagator path-integral (QUAPI) method 



are from M.Thorwart's article. 
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FIG. 3: The dynamics of non-equilibrium correlation function P{t). The parameters are: a = 0.01, 
r = 0.02. 
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FIG. 4: The dynamics of non-equilibrium correlation function P(t) for the on-resonance case. The 
parameters are: A = F = 0.02. 
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U2 corresponding to A and f2. First figure shows ooi is the dominant frequency for the off-resonance 
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FIG. 6: The validation of Shiba's relation. 
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FIG. 8: The dynamics of non-equilibrium correlation function P{t) for off-resonance case, main 
plot: P{t) is very like a over damping curve. The parameter is a = ac = 0.49997. Inset: Another 
P{t) for comparison, where one inherent frequency is dominating the dynamics of non-equilibrium 
correlation function. The parameter is: a = 0.01. Other parameters are: A/Cl = 0.5, T = 0.02. 
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FIG. 9: The dynamics of non-equilibrium correlation function P{t) for on-resonance case, main 
plot: P{t) is very like a over damping curve. The parameter is a = Oc = 0.156. Inset: Another 
P{t) for comparison, where two inherent frequency are dominating the dynamics of non-equilibrium 
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